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Abstract 

We reconsider the toy model studied in [1] of a spatially closed Friedmann-Lemaitre 
universe, driven by a massive scalar field, which deflates quasi-exponentially, bounces and 
then enters a period of standard inflation. We find that the equations for the matter 
density perturbations can be solved analytically, at least at lowest order in some "slow-roll" 
parameter. We can therefore give, in that limit, the explicit spectrum of the post-bounce 
perturbations in terms of their pre-bounce initial spectrum. Our result is twofold. If the 
pre-bounce growing and decaying modes are of the same order of magnitude at the bounce, 
then the spectrum of the pre-bounce growing modes is carried over to the post-bounce 
decaying modes ("mode inversion"). On the other hand, if, more likely, the pre-bounce 
growing modes dominate, then resolution at next order indicates that their spectrum is 
nicely carried over, with reduced amplitude, to the post-bounce growing modes. 
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I. Introduction 



Four dimensional, general relativistic and bouncing Friedmann-Lemaitre models, to 
which string-inspired "pre-Big Bang" [2] and "ekpyrotic" or "cyclic" [3] universes could 
reduce within some effective theory limit, have recently attracted renewed interest. A still 
debated issue however is how the spectrum of initial, pre-bounce, matter fluctuations is 
modified by the bounce (see e.g. [4]). 

As a warm up exercise, a toy model of a bouncing universe was studied in [1] (see 
also [5]) : a spatially closed Friedmann-Lemaitre model, driven by a massive scalar field, 
which deflates quasi-exponentially, bounces and then enters a period of standard inflation. 
Unfortunately, no definite prediction on the post-bounce spectrum of perturbations was 
reached, the main reason being the singular behaviour of the evolution equation in the 
bounce region. 

In this Note, we reconsider this simple model and rewrite the evolution equation for 
the matter perturbations in a well-behaved form. Having done so, we are able to solve 
it analytically, at least at lowest order in some "slow-roll" parameter, that is, when the 
pre and post-bounce quasi-de Sitter periods are long enough. We shall hence obtain the 
explicit spectrum of the perturbations when they exit the Hubble radius during post- 
bounce inflation, in terms of their initial spectrum when they enter the Hubble radius 
during pre-bounce deflation. 

As we shall see, two cases will arise. If the pre-bounce growing and decaying modes are 
of the same order of magnitude at the bounce, then "mode mixing" will turn out to reduce 
to mode inversion : the pre-bounce spectrum of the cosmologically interesting growing 
modes is carried over to the post-bounce decaying modes, and hence soon lost. (As for 
the post-bounce surviving growing modes, they inherit the pre- "big-bang" decaying mode 
spectrum, which, usually, is unfortunately blue, see [1-3].) 

On the other hand, if the pre-bounce decaying modes have become negligible at the 
bounce, then resolution at next order of the perturbation equation indicates that the 
pre-bounce growing mode spectrum (usually scale invariant) is nicely carried over to the 
surviving, cosmologically interesting, post-bounce growing modes with no modification, 
apart from an overall reduction factor. 

We shall conclude by a few remarks on the genericity of the result and the validity of 
the slow- roll approximation which was made to yield it. 



Consider a spatially closed, homogeneous and isotropic universe with line element : 
ds 2 = —dt 2 + S 2 (t)dQ^ where t is cosmic time, S(t) the scale factor and dVt 2 the line 
element of a 3-dimensional unit sphere. If matter is just a scalar field <p with mass m 
the Einstein equations reduce to the homogeneous Klein-Gordon and Friedmann-Lemaitre 
equations : 



II. The background 





(2.1) 



where a dot denotes differentiation with respect to cosmic time, where H = S / S is the 
Hubble parameter and where k is Einstein's constant. 
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The system of equations (2.1) has been thoroughly studied, in particular in [6-8]. We 
shall retain here that there exist ranges of initial conditions for which the scale factor has 
a minimum. We shall restrict ourselves to the case when such a bounce occurs, at t = 
without loss of generality, and set the initial conditions there : <p(0) = tpo and </?(0) = (po 
(the initial condition for S follows from the fact that S(0) = 0). Introducing rescaled 
initial condition, time, scalar field and scale factor as 

00 = \ tVo , T = 4> Q mt , 4>= — , a = 4> mS, (2.2) 
V 6 (p 

as well as the auxiliary function z(r) = — a/¥ <&, the system (2.1) becomes : 



d(j) z dz 3zda da / / z 2 \ 



where the plus sign holds after the bounce (r > 0) and the minus sign before. As for the 
initial conditions, they become 

0(0) = 1 , z(0) = z , o(0)= j_= (2.4) 

v 1 + z om 

(with zq = ^ ^J~J^ cpo). As one can see from (2.3), the solution for r < corresponding to the 
set of initial conditions (0o, zq) can be obtained from the solution for r > corresponding 
to the set (0o, —zq) by means of the transformation 

a(r, zq) = a(-r, -z ) , 0(r, z ) = 0(-r, -zo) , ^(t, ^o) = -^(-t, -^o) 4>o nx ed . (2.5) 

Now, if the standard 1 conditions for post-bounce inflation are imposed, that is if 

00 > 1 and |z o |<0o (2.6) 

then the solution of the system (2.3) can be approximated, at zeroth order in the "slow-roll" 
parameter 1/0q, by (see [1] who limited themselves to the case zq = 0) : 

i , zq 1 sinhr , . 2 rtS , n _ s 

a^coshr , 0~1 , 2~ ^ o— (cosh r + 2) . (2.7) 

cosh r 3 cosh r 

By comparison with direct numerical integration of (2.3), one sees that (2.7) is a good 
approximation to the exact solution for and a as long as |r| <C 0ch an d a good approx- 
imation for z on the much wider range |r| 0q, that is as long as the universe is well 
within the two, pre and post-bounce, dust-like eras. 2 



1 See [7] for fine-tuned values of 0o of order 1 which also yield inflation. See [8] for 
a proof that all models recollapse. Generic values of 0o of order 1 yield small universes 
which soon recollapse. Hence the toy model considered here is inappropiate to describe a 
post-bounce universe which does not inflate and immediately enters a radiation era. 

2 The analytical solution at next order in 1/0 2 , can easily be obtained, see [1], but will 
not be used in this Note. 
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III. The evolution equation for the scalar perturbations 

We consider now the perturbed metric : ds 2 = -(1 + 2$)dt 2 + S 2 (t)(l - 2^)dQ§ and 
the perturbed scalar field if(t) + x- In Fourier space, the "scalar" perturbations <E> n , \l/ n 
and Xn are functions of time and of the Eigenvalues n of the Laplacian on the 3-sphere 
(defined as A/ n = — n(n + 2)/ n , n E N and n > 2). The (kl), (Ok) and (00) components 
of the linearized Einstein equations then are, respectively (see, e.g. [9]) : 



+ 



i> n + H<S> n =^-Xn 

3if$ n k ( dV 



(3.1) 



where k 2 = n(n + 2) — 3K, with K = 1, and where, in the toy model we consider here : 
F(<^) = \m 2 ip 2 . As is well-known [9], the last equation can be rewritten, using the two 
constraints as : 



<!>„+ \ lH+\ d f\ $ n + 

ip d(p ) 



k< 



-5K f HdV 

-s 2 - + 2 { KV + Jd^ 



®r, = 



(3.2) 



Another form of that equation is easily found to be, see [9] : 



+ [k 2 - W(rj)] u n = with 



j ii 



,12 



U 



, -^ n and W( V ) = -^ t + 2^ + k^- (3.3) 



if 



if 



where a prime denotes differentiation with respect to conformal time r] — related to cosmic 
time t by Sdrj = dt. 

It is clear that none of the forms (3.1) (3.2) or (3.3) is suitable for integration when 
<f goes through zero (which is necessarily the case if there is to be quasi-de Sitter regimes 
before and after the bounce, see e.g. equation (2.7) for z oc if). 3 Now, it is easy to 
check that, at least when V(<p) = \m 2 if 2 ^ they can be put into the strictly equivalent, 
well-behaved form 



d(a$ r 



dr 
dg n = 
dr a 5 (f> 



— A 

r 

H To" ( 3 — 



k' 2 



a 2 (j) 2 



with A r 



k 2 $ 7 



(3.4) 



where a(r), 0(r) and z(r) solve the background equations (2.3). (The first equation is 
nothing but the (Ok) linearized equation and the second is a rewriting of the (00) one in 



3 The authors of ref [1] solved the perturbation equation (3.2) (in the particular case 
zq = 0) using delicate numerical matching techniques. 

4 The result can easily be extended to any potential of the form V(f) = cq + c\if n . 
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terms of the suitably chosen auxiliary function g n . 5 ) Once (3.4) is solved and $ n and g n 
known, then the other scalar perturbations are given by 

= and ^YTo = ^ = { Zgn ~^J • (3 - 5) 

Note that, in view of the symmetry properties of the background solution, see (2.5), the 
solution of (3.4) is such that 

$ n (r, zq, $ n (0), g n (0)) = $ n (-r, -zq, $ n (0), -g n (0)) . (3.6) 



IV. Relating the pre and post-bounce spectra at lowest order 
in i /on 

When the conditions (2.6) on the initial conditions are met, the term in 1/4>q in 
the system (3.4) can, at lowest order, be ignored. The evolution equation for the scalar 
perturbations thus simplifies into 



d(a$>, 
dr 



k " A 
~ A 

S± r 



Gr 



dg r 

dr a 5 



with, recall : A n = a 3 ( zg n — 



a 2 (f) 



(4.1) 



where, at the same approximation, the background functions a(r), 4>(r) and z(t) are given 
by (2.7). Differentiating A n once, one finds : 



dA n 
dr 



-5r n cosh r. 



(4.2) 



Differentiating again, one gets the following, closed, equation for A n (or, equivalently for 

Xn) ■ 



d 2 A r 
~d^ 



dA n 



— 3tanhr— 1 o~A n ~ . 

dr cosh t 



(4.3) 



Recalling that k 2 = (n — l)(n + 3), the general solution of equation (4.3) is a sum of 
even and odd functions : 

A n = a n A^+p n A^ (4.4) 
where (a n , f3 n ) are constants of integration, and 
• for t > : 



A { y = cosh 3 rF 



n + 2 n _1 
2 ' 2 ' 11 



1 



cosh r 



= F 



l—n n+3 5 
2 ' 2 '25 



1 



cosh T 



(4.5) 



5 In the late-time dust-like era when ip and <p oc (/» go periodically through zero, but H 
remains positive, another well-behaved form is required and was given in [10]. 
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for t < : 
Ai 1] = ± cosh 3 r F 



n+2 n _1 
2 ' 2 ' 21 



1 



cosh T 



, A^= T F 



1 — 71 n.+3 5 
2 ' 2 '2' 



1 



cosh T 



(4.6) 



where the upper vs lower signs hold for n even n odd, and where F[o,6,c, x] is the 
hypergeometric function (usually denoted 2Fl[a,b,c, x]). 

The function A n being known, the scalar perturbation $ n follows from (4.1-2) and 
the approximate background solution (2.7). It reads : 



$(!> 2 ) 



1 



1 



(n — 1) (n + 3) cosh r 



with 



cL4 



(1.2) 



(4.7) 



(n- l)(ra + 3) 



v^M 1,2) + 



2cL4 



(1.2) 



3 e&r 



x vl — a; Vl — #(1 + 2x) =F 3zq£ 



3/2 



where x = a 2 ~ cos | i 2 - and where the upper vs lower sign hold for (r > 0) (r < 0). At 
the bounce : 



r=0 — 



,(2) 



(n- l)(ra + 3) 



(n + 1) cos — + ^(^ + 2) sin — 



n r=0 



(ra-l)(ra + 3) 



sin ^ z n7r 

— — r COS 

n(n + 2) n + 1 2 



(i) . , . nn (2) 

Q I n = -n fi + 2 sin — , q\ n 

^>n t=0 v 1 / O ' r=0 



3 nn 

cos — — 

n + 1 2 



(4.8) 



(4.9) 



Writing down the explicit expression of $ n in terms of hypergeometric functions is not 
particularly illuminating : suffice it to say that it is a good approximation of the numerical 
solution of the exact equations (2.3) (3.4) in the range |r| <C (j>o, and that it tends to 
constants for large |r| (in practice |r| bigger than a few unities) and oscillates in the 
bouncing region more and more as n grows bigger. 

We are now in a position to relate the pre and post-bounce spectra, at the order 
considered, that is the lowest in the slow-roll parameter 1/4>q. 

Using the following asymptotic expansions of hypergeometric functions, 

F[^,^,l,x] = l + 0(x) , F x] = l + n(n o + 2) x + 0(x 2 ) , (4.10) 



the asymptotic behaviours of $ n (in practice for |r| bigger than a few unities) are readily 
obtained from (4.7) (recalling that x ~ a~ 2 ) : 



D n 

• post-bounce region : $ n ~ G™ t H ^— 1 with 



a 



/-in 
'-''post 



nn _ 
-^post — 



3 

(n-l)(n + 3) 



(4.11) 



Qn -^pre i g 

• pre-bounce region : $ n ~ e + -J^- with ^ /j n + 3^ a n ( 4 ' 12 ) 

I pre " (n-l)(n + 3) 

where the upper vs lower signs hold for even vs odd n. 

At lowest order in the slow- roll parameter 1 the pre and post-bounce spectra are 
thus very simply related : 

t::Hfl ")Ki) +o(i/ «- <4 - i3) 

In the toy model and at the approximation considered here, "mode-mixing" therefore 
reduces to a simple "mode inversion", that is the pre-bounce spectrum of the growing 
modes is carried over to the post-bounce decaying modes (and vice- versa). In other words : 
the pre-bounce spectrum of the cosmo logically interesting growing modes is carried over to 
the post-bounce decaying - modes, and hence soon lost. On the other hand, the post-bounce 
surviving growing modes, inherit the pre-bounce decaying mode spectrum, which, usually, 
is unfortunately blue [1-3]. 

There is however a case when the result does not hold, to wit a n = 0, that is when the 
pre-bounce decaying modes have become insignificant when reaching the bouncing region. 
This is a case of physical interest, to which we now turn. 

V. Validity of the approximation and next order in 

Let us first look at the structure of the "tranfer matrix" at next order in the slow-roll 
parameter. Including 1/4>q corrections will yield 



GUA „ ( c(n)/$ ±1 \ [G- re 



£>pW~V =Fl d(n)/cf>lj \D- Ie 



(5.1) 



and will not change significantly the zeroth order result, if -Dp re and G™ ve are of the same 
order of magnitude. If, now, the initial conditions on D™ re and G™ Te are such that 

which, in view of (4.12), is the case when a n is vanishingly small so that the pre-bounce 
decaying modes have become insignificant when reaching the bouncing region, then the 
pre and post-bounce spectra become related by 

^post — ,2 ^pre > -^post — T^pre • \°-°) 

This is the physically relevant case studied numerically in [1] where Gp OSt is found to be 
non zero even when D™ Ye = 0. The authors of [1] did not however give the n-dependance 
of the constant c(n) and hence left open the question of how G™ ost was related to G™ re . 
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The constant c(n) can however be estimated as follows. 

The exact equations for the perturbations are (3.4) where the background functions 
solve (2.3). In the previous Section we solved them at zeroth order in 1/</>q, that is we 
ignored the 1/4% term in (3.4) and used for the background functions the zeroth order 
approximation (2.7). To consistently iterate them at next order one should : 

1. keep the 1/4% term in (3.4), replacing <£> n by the zeroth order solution ; 

2. use for the background functions the first order approximation of (2.3). 

To estimate c(n), we shall however ignore step 2, for the following reason : treating 
the background at first order in the slow-roll parameter introduces logarithmic corrections 
in the solutions of the perturbation equation which render more difficult the numerical 
extraction of the post-bounce growing mode. Of course, the difficulty can be overcome, 7 
but is perhaps not worth the effort as it seems a reasonable guess that the inclusion of the 
1/00 correction for the backgound should not affect the n-dependence of the perturbation 
spectrum. 

We therefore integrated numerically the set of equations (3.4) where, in the 1/4% term, 

(2) 

we replaced $ n by the pre-bounce purely growing zeroth order solution, that is <3>n , see 
(4.5-7), but where we used for the background functions the zeroth order approximation 

(2-7). ' , , , , 

We chose the initial conditions at the bounce, given by (4.8-9) with a n = 0, (3 n = 1. 
With those initial conditions, integration yields a solution <£>^ er which, when 4o — * oo, is 
nothing but the analytical solution obtained in the previous Section, that is a mode 
which is exponentially growing before the bounce and exponentially decaying after the 
bounce. For large but finite value for 4% on the other hand, $^ er no longer vanishes in the 
asymptotic regions but tends to small constants, (Cp OSt , Cp re ), which scale, as they should, 
as 1/4%- Let us therefore introduce the rescaled constants (Cp OSt = 4o c p OS t' Cpre = 0o c pre) 
which are independent of 4%- 

Consider now the linear combination 

C n 

^shooting = $iter ± 3-gp$^ (5.4) 

where the upper vs lower signs hold for n even vs n odd. Since — ► Tl/3 in the 
pre-bounce asymptotic region (see (4.12)), we have 

• pre-bounce region : $^ hootm s _ > (5.5) 

so that $^ hootin g i s the purely pre-bounce growing mode at the approximation considered. 
In the post-bounce region on the other hand (see (4.11)) : 

(jn -j- Qn 

• post-bounce region : $f ootin g _> post J 2 pre (5.6) 

4>Q 

7 Indeed, the task of integrating numerically the background and the well-behaved per- 
turbation equations (3.4) is in principle straightforward. 
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Therefore, from (4.12) and (5.3) 



G£ ost ^ T n G£ re with T n = ±^(C; ostT C; re ). (5.7) 

Numerical integration gives, pour each n (and each value of the parameter zq entering 
the background zeroth order solution (2.7)) the values of the constants Cp OSt and Cp re and 
it turns out, remarkably, that, for large n, T n does not depend on n. Hence the pre-bounce 
spectrum encoded in Gp re is nicely carried over to the post-bounce growing mode Gp OSt , 
although with an amplitude reduced by the overall factor 1/4>q. 

VI. Concluding remarks 

In the very simple toy model of a bouncing universe that we studied in this paper, and 
at lowest order in the slow-roll parameter, the spectrum of the pre-bounce perturbations is 
carried over through the bounce with a simple (although unfortunate) inversion of modes. 
That analysis however breaks down if the pre-bounce decaying modes are negligible in the 
pre-bounce region. In that case the analysis must be pushed to next order in the slow-roll 
parameter with the neat indication that the large n post-bounce growing modes inherit 
without distortion the pre-bounce growing mode spectrum. 

These results do not depend on the value of <fi at the bounce. (Indeed numerical 
integration indicates that the value of zq only affects the overall amplitude of the tranfer 
factor T n .) 

It would be surprising if they depended crucially on the particular potential (V(ip) = 
^m 2 ip 2 ) chosen for the scalar field, as long as there exist quasi-de Sitter regimes before 
and after the bounce. They should not be spoilt either when treating more carefully the 
logarithmic corrections, but that point certainly deserves further attention. 

Of course, it is not clear whether current string-inspired bouncing universes ("pre- 
Big-Bang", "ekpyrotic" or "cyclic") can fit such a simple framework. If not, toy models as 
the one considered here become irrelevant. 



Acknowledgements. N.D. is grateful to the organisors of the stimulating cosmology 
workshop which took place in Cambridge in July 2003, and thanks in particular Neil 
Turok and Paul Steinhardt for very interesting discussions. 

References 

[1] C. Gordon, N. Turok, Phys. Rev. D 67 (2003) 123508 

[2] G. Veneziano, Phys. Lett. B 265 (1991) 287; M. Gasperini, G. Veneziano, Phys. Rep. 
373 (2003) (and references therein) 

[3] J. Khoury, B.A. Ovrut, P.J. Steinhardt, N. Turok, Phys. Rev. D 64 (2001) 123522; J. 
Khoury, B.A. Ovrut, N. Seiberg, P.J. Steinhardt, N. Turok, Phys. Rev. D 65 (2002) 
086007; N. Turok, P.J. Steinhardt, Phys. Rev. D 65 (2002) 126003; J. Khoury, B.A. 



9 



Ovrut, P.J. Steinhardt, N. Turok, Phys. Rev. D 66 (2002) 0406005 (and references 
therein) 

[4] R. Brustein, M. Gasperini, M. Giovaninni, V. Mukhanov, G. Veneziano, Phys. Rev. 
D 51 (1995) 6744 ; M. Gasperini, N. Giovaninni, G. Veneziano, Phys. Lett. B 
259 (2003) 113; Brandenberger, F. Finelli, JHEP 0111 (2001) 056; R. F. Finelli, 
R. Brandenberger, Phys. Rev. D 65 (2002) 103522; A.J. Tolley, N. Turok, |Ke]> 



th/0204091| ; D.H. Lyth, Phys. Lett. B 524 (2002) 1; J.C. Hwang, Phys. Rev. D 65 



(2002) 063514; S. Tsujikawa, Phys. Lett. B 526 (2002) 179; P. Peter, N. Pinto-Neto, 
Phys. Rev D 65 (2001) 123513 and Phys. Rev D 66 (2002) 063509; J. Martin, P. 
Peter, N. Pinto-Neto, D.J. Schwarz |hep-th/0112128| , 0204222, 0204227; J. C Fabris, 
R.G. Furtado, P. Peter, N. Pinto-Neto, Phys. Rev D 67 (2003) 124003, P. Peter, N. 
Pinto-Neto, |hep-th/0306005| ; J. Martin, P. Peter, Phys. Rev D 68 (2003) 103517; R. 
Durrer, F. Vernizzi, hep-th/0203375| ; L. Allen and D. Wands, |astro-ph/ 0404441 
[5] J. Hwang, H. Noh, Phys. Rev. D 65 (2002) 124010 

[6] H. Nariai, Prog. Theor. Phys. 46 (1971) 433; H. Nariai, K. Tomita, Prog. Theor. 
Phys. 46 (1971) 776; L. Parker, S.A. Fulling, Phys. Rev. D 7 (1973) 2357; H. 
Nariai, Prog. Theor. Phys. 51 (1974) 613; A.A. Starobinski, Pisma A.J. 4 (1978) 
155; V.N. Melnikov, S,V. Orlov, Phys. Lett. A 70 (1979) 263; J.D. Barrow, R.A, 
Matzner, Phys. Rev. D 21 (1980) 336; J. B. Hartle, S.W. Hawking, Phys. Rev. D 
28 (1983) 2960; S.W. Hawking, J.C. Luttrell, Nucl. Phys. B 247 (1984) 250; D.N. 
Page, Class, and Quant. Grav. 1 (1984) 417; V.A. Belinsky, L.P. Grishchuk, Ya.B. 
Zel'dovich, I.M. Khalatnikov, J. Exp. Theor. Phys. 89 (1985) 346; M. Visser, Phys. 
Lett. B 349 (1995) 443; A. Yu Kamenshchik, I.M. Khalatnikov, A. V. Toporensky, 
Int. Journ. Mod. Phys. D 6 (1997) 673; Int. Journ. Mod. Phys. D 7 (1998) 129; 
I.M. Khalatnikov, A. Yu Kamenshchik, Phys. Rep. 288 (1997) 513; C. Molina-Paris, 
M. Visser, Phys. Lett. B 455 (1999) 90; S. Gratton, A. Lewis, N. Turok, Phys. Rev 
D 65 (2002) 043513 

[7] N.J. Cornish, E.P.S. Shellard, Phys. Rev. Lett. 81 (1998) 3571 

[8] H.J. Schmidt, Astron. Nachr. 311 (1990) 99 

[9] J.M. Bardeen, Phys. Rev. D 22 (1980) 1882; M. Sasaki, Prog. Thor. Phys. 70 (1983) 
394; V.F. Mukhanov, H.A. Feldman, R.H. Brandenberger, Phys. Rep. 215 (1992) 
203; J.M. Stewart, Class. Quant. Grav. 7 (1990) 1169 
[10] N. Deruelle, C. Gundlach, D. Polarski, Class. Quant. Grav, 9 (1992) 137 



10 



